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Abstract. We perform a detailed analysis of first order Sobolev-regular infinitesimal isometries 
O-l on developable surfaces without affine regions. We prove that given enough regularity of the 

surface, any first order infinitesimal isometry can be matched to an infinitesimal isometry of an 
arbitrarily high order. We discuss the implications of this result for the elasticity of thin devel- 
00 . opable shells. 
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?5 ■ 1. Introduction 

The derivation of asymptotic theories for thin elastic films has been a longstanding problem in 
the mathematical theory of elasticity [2]. Recently, various lower dimensional theories have been 
>^ • rigorously derived from the nonlinear 3 dimensional model, through F-convergence [5j methods. 

I Consequently, what seemed to be competing and contradictory theories for elastica (rods, plates, 

■ shells, etc) are now revealed to be each valid in their own specific range of parameters such as 

material elastic constants, boundary conditions and force magnitudes [211 [HI [U [281 H] • Iii this line, 
Priesecke, James and Miiller gave a detailed description of the so called hierarchy of plate theories 
in [9], corresponding to distinct energy scaling laws in terms of the plate thickness. Similar results 
! have been established for elastic shells [22\ [71 [23l [Mj [25] , however the description is still far from 

being complete. 

In [26j Lewicka and Pakzad put forward a conjecture regarding existence of infinitely many 
small slope shell theories each valid for a corresponding range of energy scalings. This conjecture, 
based on formal asymptotic expansions, is in accordance with all therigorously obtained results for 
5^ , plates and shells. It predicts the form of the 2 dimensional limit energy functional, and identifies 

the space of admissible deformations as infinitesimal isometries of a given integer order > 
determined by the magnitude of the elastic energy. Hence, the influence of shell's geometry on 
its qualitative response to an external force, i.e. the shell's rigidity, is reflected in a hierarchy of 
functional spaces of isometries (and inflnitesimal isometries) arising as constraints of the derived 
theories. 

In certain cases, a given A^th order infinitesimal isometry can be modified by higher order 
corrections to yield an infinitesimal isometry of order M > N, a property to which we refer to by 
matching property of infinitesimal isometries. This feature, combined with certain density results 
for spaces of isometries, causes the theories corresponding to orders of infinitesimal isometries 
between N and M to collapse all into one and the same theory. Examples of such behavior are 
observed for plates |9], where any second order infinitesimal isometry can be matched to an exact 
isometry (and hence to an Mth-order isometry for all M G N), and for convex shells [25j, where 
any first order infinitesimal isometry enjoys the same property. The effects of these observations 
on the elasticity of thin films are drastic. A plate possesses three types of small-slope theories: the 
linear theory, the von Karman theory and the linearized Kirchhoff theory [9j, whereas the only 
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small slope theory for a convex shell is the linear theory [25]: a convex shell transitions directly 
from the linear regime to fully nonlinear bending if the applied forces are adequately increased. 
In other words, while the von Karman theory describes buckling of thin plates, the equivalent 
variationally correct theory for elliptic shells is the purely nonlinear bending. 

In this paper, we focus on developable surfaces (without flat regions). This class includes smooth 
cylindrical shells which are ubiquitous in nature and technology over a range of length scales. An 
example of a recently discovered structure is carbon nanotubes, i.e. molecular-scale tubes of 
graphitic carbon with outstanding rigidity properties [TT]: they are among the stiffest materials 
in terms of the tensile strength and elastic modulus, but they easily buckle under compressive, 
torsional or bending stress [18j. The common approach in studying buckling phenomenon for 
cylindrical tubes has been to use the von Karman-Donnel equations [27^ I13j . However, as we 
establish here, the proper sublinear theory for this purpose is, again, the purely nonlinear bending 
theory. It seems likely that the von Karman-Donnel equations could be rigorously derived and 
valid in another scaling limit, e.g. when the radius of the cylinder is very large as the thickness 
vanishes. 

The key ingredient of our analysis is a study of VF^'^ first order isometrics on developable 
surfaces, of which we address regularity, rigidity, density and matching properties. Our results 
depend on the regularity of the surface and a certain mild convexity property. More precisely, we 
establish that any c2Af-i,i j-gg^ia,r first order infinitesimal isometry on a developable C^^'^ surface 
with a positive lower bound on the mean curvature, can be matched to an A^th-order infinitesimal 
isometry. Combined with a density result for VF^'^ first order isometrics on such surfaces, we prove 
that the limit theories for the energy scalings of the order lower than collapse all into the 

linear theory. Our method is to inductively solve the linearized metric equation symVw = B on 
the surface with suitably chosen right hand sides, a process during which we lose regularity: only 
if the surface is C°° we can establish the total collapse of all small slope The importance of the 
solvability of symVw = B has been noted in [33], see also [10] and [1] in regards to relations 
of rigidity and elasticity. We remark here that a surface, e.g. a regular cylinder, may well be 
ill-inhibited according to the definition by Sanchez-Palencia in [33] and yet satisfy the adequate 
matching properties resulting in the aforementioned collapse. 

Our analysis can be generalized to piecewise smooth surfaces which satisfy the convexity prop- 
erty as above, on each component. The question of existence of a developable or non-developable 
surface with no flat regions which shows a different elastic behavior (i.e. validity of an intermediate 
theory between the linear one and the purely nonlinear bending) remains open. 

Finally, a word on the developability property of surfaces of vanishing Gaussian curvature, 
from which the term developable is derived, is to the point. For each point on such a surface there 
exists a straight segment passing through it and lying on the surface, in a characteristic direction. 
Developable surfaces are also locally identified with isometric images of domains in M?; this last 
property and, in particular, the developable structure of isometrics of flat domains is heavily 
exploited in this paper. Such structure was established for isometrics in [T2], for isometrics 
with total zero curvature in [3J4 Chapter II], ^32^ Chapter IX] and for W'^''^ isometrics in [T9l [30]. 
In |30j Pakzad proved that any M^^'^ isometry on a convex domain can be approximated in strong 
norm by smooth isometrics, and in [29] the boundary regularity was discussed. Lately, Hornung 
systematically represented and generalized these results in jl4[ [T5] [T6] whose terminology we will 
adapt for the sake of simplicity and completeness. The above mentioned results have had other 
applications in nonlinear elasticity [31 [20l [T7] . 

The paper is organized as follows. In section [21 we introduce and review preliminary facts about 
developable surfaces. In section [3l we study the linearized equation of isometric immersion of a 
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developable surface and prove existence of a solution operator with suitable bounds in Theorem 
13.11 We proceed to study the space of W'^''^ first order infinitesimal isometrics and prove a 
compensated regularity and rigidity property of such mappings in Theorem 14.11 of section [H 
We then use these results to prove in Theorem 15.21 that given enough regularity of the surface, 
any first order infinitesimal isometry can be matched to a higher order one. Combined with a 
straightforward density result in Theorem 15.31 we are finally able to derive the main application 
of this paper in elasticity theory, namely Theorem 16.21 which is the counterpart to Theorem 16.11 
for deducing the F-limit of 3 dimensional nonlinear elasticity for thin developable shells. 

Acknowledgments. P.H. was supported by EPSRC grant EP/F048769/1. M.L. was partially 
supported by the NSF grants DMS-0707275 and DMS-0846996, and by the Pohsh MN grant 
N N201 547438. M.R.P. was partially supported by the University of Pittsburgh grant CRDF- 
9003034 and by the NSF grant DMS-0907844. 

2. DEVELOPABLE SURFACES 

Let C be a bounded Lipschitz domain, and let u G be an isometric immersion 

ofn into R^: 

diu{x) ■ dju{x) = 5ij for all x £ ^l. 

A classical result jl2j asserts that, away from affine regions, such u must be developable: the 
domain il. can be decomposed (up to a controlled remainder) into finitely many subdomains on 
which u is afhne and finitely many subdomains on which u admits a line of curvature parametriza- 
tion, see e.g. Theorem 4 in [15]. In light of this result, it is natural to restrict ourselves to the 
situation when can be covered by a single line of curvature chart. 

We now make the above more precise. Let T > 0, let F G C^''^{[0,T],R^) be an arclength 
parametrized curve, and let G C'^'^([0, T]) be positive functions (a rationale for assuming 
Lipschitz continuity here can be found in Lemma 2.2 in [17]). Following the notation in |16) . we 
set: 

iV = (F')^ G C°'H[0,r],M2) and k = F" • iV G L°°(0, T) 
the unit normal and the curvature of F. We introduce the bounded domain: 

M,± = {{s,t) : t G (0,T), s G {- 3- (t) , s+ (t))} , 
the mapping $ : Mg± — )• M? given by: 

<i>is,t) = r{t) + sN{t) 

and the open line segments: 

[rit)] = {rit) + sN{t):se{-s-{t),s+m- 

From now on we assume that: 

(2.1) [F(ti)] n [F(t2)] = for all unequal ti,t2 G [0, T]. 

This condition will be needed to define a developable isometry u, which maps segments [F(t)] to 
segments in M^. 

The next lemma is a consequence of Proposition 2 and Proposition 1 in [16j. 
Lemma 2.1. Given k„ G L'^{0,T), let r = {j',v,n)^ G W'^'^{{0,T), SO{^)) satisfy the ODE: 

K Kr, 

(2.2) / = I -K I r 

-Kr, 
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with initial value r(0) = Id, where we set 7(t) = Jq 7'- Define the mapping: 

(2.3) n:$(A^,±)^M^ u{<^{s,t)) = 'y{t) + sv{t) y{s,t) £ M,±. 
Then u is well defined and: 

is an isometric immersion. Moreover, Vu(4>) e C°(M^± , M^^^) with: 

(2.4) Vums,t)) =j'{t)^T\t) + v{t)^N{t) V(s,t)GM,±. 
and the functions: 

a-ij = {diu X d2u) ■ didjU 

satisfy: 

(2.5) didjU = aij{diu x d2u) Vi,j = 1,2, 

K (t) 

(2.6) aij{^{s,t)) = - — ':^^^-^r'i{t)T'j{t) for almost every {s,t) e M,±. 

1 SKyt J 



(2.7) / I / ,^M}^ds]dt < 00 



If, in addition: 

SK{t) 

then u e Ty2'2($(M^±)),M3) wii/i.- 



(2.8) /" |V2n(x)|Mx= [ 

J<I>(M+) Jo 



ds]dt. 



'*(M^±) Jo ^Js-(t) l-SK{t) 

In order to keep the hypotheses short, we make the following definition: 

Definition 2.2. A surface S CM.^ is said to be developable of class C'^'^ if there are F, N, s^, 7, 
V, n, K, Kn, ^ and u as in Iemma \2.1\ such that 



(2.9) $(M,±) = 0, 

u £ C'''^{Tl,R^) and S = u{n). 



Remark 2.3. (1) The curve 7 is a line of curvature of surface S. Hence the mapping (s, t) 1— )• 
^{t) + sv{t) is a line of curvature parametrization of S and the condition (j2.9p is the precise 
formulation of our assertion that O is covered by a single line of curvature chart. 

(2) The moving frame r is the Darboux frame on the surface S along 7. Therefore, ()2.2p 
indicates that the geodesic curvature of 7 coincides with the curvature k of its preimage 
r, and that its geodesic torsion vanishes. This is naturally expected as u is an isometry. 
In the same vein, Kn is the normal curvature of 7 on 5. 

(3) Of course, n is the unit normal vector to S and — [ojj] is the second fundamental form 
(expressed in u-coordinates). Equation (12. 6p shows that it has rank one or zero, hence the 
Gauss curvature det[ajj] is zero. 

(4) Condition (j2.7p implies that k„ = almost everywhere on the set: 



lo = {t G (0,r) : K{t) G {l/s-{t), l/s+{t)]]. 
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Now, by Proposition 10 in pL5j and in view of (12. ip we obtain: 

1 1 



(2.10) K{t) G 



for almost every t G (0,T). 



Moreover, the Lipschitz map ^ with: 

(2.11) det V^>(s,t) = -(1 - SK(t)) for almost every (s,t) G Ms±, 

is a homeomorphism from Mg± onto 0. We will frequently make the extra assumption that mean 
curvature of u be bounded away from zero: 

(2.12) |trace [aij](a;)| > for all x gQ. 
Then ()2.6p and (|2.7p imply that there is 5 > such that 



(2.13) n{t) G 



. s-{tys+{t) . 



for almost every t G (0, T). 



By Proposition 10 (iii) in [15j and the bounds (|2.13p . <1> ^ is Lipschitz as well (this assertion is 
generally false if (I2.13P is violated). 

Lemma 2.4. Assume that S is developable of class C^'^ for some k > 2 and that (j2.12p holds. 
Then k, Kn £ C^~^'^ and ^>~^ are C^~^'^ up to the boundary of their respective domains Mg± 
and rj. 

Proof. The hypothesis on S imply that aij G C'^~^'^(r2). By continuity of aij and $ and by (j2.12p . 
we may assume without loss of generality that Trace [ajj]($) > c > on M^±. In view of (|2.6p 
we have: 

(2.14) Trace [a^,]{^) = 

Since <^ is bilipschitz, this implies that the right-hand side of (12.141) is Lipschitz. As k„ > c > 
it follows that Lipschitz. Thus P', iV belong to C^'\ hence so does By (f2TT]l . (|2T3]) . 

the Jacobian of <^ is uniformly bounded away from zero on M^i, and so 4>~^ belongs to C^'^, too. 
If A; > 3 then we return to (j2.14p . apply Lemma 12.51 and argue as before to conclude that k, k„ 
are in C^'^. The conclusion follows by iteration. ■ 

In the proof of Lemma 12.41 we used the following chain rule, which is a particular case of 
Theorem 2.2.2 from [34J. 

Lemma 2.5. Let Ui,U2 C M" be two open, bounded sets and let ^ : Ui ^ U2 be a bilipschitz 
homeomorphism. Then f G W^''^{U2) if and only if f o ^ ^ W^''^{Ui). If this is the case, then the 
chain rule applies: 

(2.15) V(/o$) = ((V/)o$)V$ a.e. in Ui. 

3. Equation symVw = i? on developable surfaces 

In this section we let S" be a developable surface of class C^'^. By n : — )• we denote the 
unit normal to S which satisfies |n(x)| = 1 and n(u(x)) = diu{x) x d2u{x) for all x G ^2, and 
^ili't)) = iT-it) for all t G (0, r). By 11 = Vn we denote the second fundamental form of S, defined 
symmetric bilinear form by: 

n(j')(r, rj) ='n- drn, \/ T,ri £ TpS, p £ S, 

where TpS is the tangent plane to S at point p, so that: 

Il{u{x)){diu{x) , dju{x)) = —aij{x) for all 3; G il. 
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We continue to assume ()2.12p . Hence: 

(3.1) |K„(t)|>c>0 VtG[0,r]. 

For a given symmetric bilinear form B £ C^'^{S),M?^'^, we want to solve a first order PDE: 

(3.2) symVw = B, 

on 5, where to : 5" — )• is a displacement field and the expression symVw in the left-hand side 
is the following bilinear form acting on the tangent space of S: 

sjmVw{p){T, ri) = ]^{drw{p) ■ rj + d^w{p) • r), V r, e TpS, p£ S. 

1. We shall write wn to denote the scalar product w ■ n, and we decompose w as follows: 

w = wtan + {wn)n. 

Hence wtan{p) G TpS for all p £ S. We define the pulled back maps ws = {wn) o u and w' = 
{wtan ° u)'^\/u, as Well as the pulled back form: 

Bij{x) = B{u{x)){^^u{x),^ju{x)) ^x £ n Vi,j = l,2. 

Using (j2.5p and recalling that wtan is tangent to S, we calculate: 



dju{x) ■ Vw{u{x))diu{x) = di{w{u{x))) ■ dju{x) = di(^w{u{x)) ■ dju{x)^ — w{u{x)) ■ dfju{x) 

= diWj - wsaijix). 

Hence (j3.2p can be written in terms of the pulled back quantities as the following matrix equality: 
(3.3) [Bij] = symVui' - wslatj], 

where now sym \/w' is understood in the usual way, with respect to the standard Euclidean 
coordinates in il. 

2. Recalling that the condition for a matrix field B to be of the form B = symVw for some 
vector field zZ; on is equivalent to curl"^curl5 = (in view of Q being simply connected), the 
equation (13. 3j) becomes: 



(3.4) 



curl curl [Bij] = —curl curl (^tDa [ojj 

= —W3 cur 1"^ curl [aij] — 2V^W3 ■ curl [uij] — cof V^t^s : [a 
Notice now that: 

dfiU ■ — 5^2^ ■ din 



curl a. 



-curl 



di2U ■ n 822^ ■ n 



di2U ■ d2n — d22U ■ din 



because di{n o u) £ TpS and d'fjU • d^u = by Lemma [27T1 Hence: 

(3.5) e = curFcurl [Bij] 
belongs to L°°. The problem ()3.4p becomes: 

(3.6) 6 = —cof V^ws : [aij] in fl. 
3. By Lemma |2. II we have: 



(cof {V'ws) : [a,,]) ms,t)) = -^^dUwms , t))) ■ 
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Consequently, problem ()3.6p is equivalent to: 

(3.7) dl{w3ms,t))) = ^ems,t)) for all {s,t) G M,±. 

The above calculations show that in order to solve (j3.2|) . it is sufficient and necessary to solve 
the ODE (13. 7p for and then recover w' from (13. 3p . Moreover, the solution (w'jW^) is unique 
after choosing the boundary conditions w^{^(t,0)) and ds{w^{^(t,0)), where uniqueness of w' is 
understood up to afHne (linearized) rotations of the form A{s,t) + 6, ^ G so{2), b G M?. 

Theorem 3.1. Assume that S is developable of class C^'^ and satisfies (I2.12p . and let a G (0, 1). 
Then there exists a constant C such that the following is true. For every symmetric bilinear form 
B G C^'"'^(S', M^^^) there exists a solution w = Wtan + {wn)n with Wtan G C^'" and (wn) G of: 

(3.8) symVw = symVwtan + {wn)Il = B 
satisfying the bounds: 

(3.9) Iktanllco." + ||(lf"^)||oo < C\\B\\ci,l. 

If, in addition, S G C*^^^'^ and B G C^^^^"^ for some k>l, then 

(3.10) llwtanllc*^! + ll('«^^)llcfe-i.i ^ C\\B\\ck+i,l. 

Proof. 1. Assume first the minimal regularity u G C^'^ and B G C^'^ , so that = curl^curl[i?jj] G 
L°°(r2,M). Solving (jST]) by integrating twice in s from t(;3(<I>(t, 0)) = and dsW^{^{t,{))) = 0, we 
obtain: 

(3.11) ll(^i'^)l|L- < c^II^IIl- < C'll^llci,!. 
Solving now (13. 3p for Wtan so that: 

skew -f V{wtan ou) =Q and -f (lytan o n) = 0, 

(where for a square matrix P, its skew-symmetric part is denoted skew P = |(P + P^)), we 
obtain by means of Korn's inequality, for any p > 1: 

\\Vwtan\\LP < C\\symV{wtanOu)\\LP < C {\\[Bij]\\L^ + || (ii;n) || L- ) < C\\B\\ci,i, 

where C may depend on p. Combining with the Poincare inequality, we get: 

(3.12) ||w^tan||vi/i'P ^ 

By Sobolev embedding, (j3.9p follows now from (j3.1ip and ()3.12p . 

2. When B G C''+^'^ and u G C''+^'\ then G C''"^'^ by Lemma El and so by (l3TP : 

(3.13) \\{wn)\\ck-i,i < C\\e\\ck-iA < C\\B\\ck+i^, 

where C may depend on S. Recalling that V^w' can be expressed as the linear combination of 
partial derivatives of symVtf', we have from (13. Sp . (13.13P that: 

II V^Wjarallc'-— 2.1 ^ C'l|-S|lc'=+i'l5 

which implies p.lOp in view of (j3.9p and (j3.13p . ■ 

Proposition 3.2. Assume that S is developable of class C^^"^'^ , satisfying 112. Also assume 

that B = sym((V(/))"^(V'(/')) where 4>,il) G C'^^"'^'^(S', M^). Then w as obtained in Theorem \3.1\ 
satisfies: 

\\wtan\\ck.i + \\wn\\ck-i.i < C||V'||cfc+i>i||</'llcfe+i.i- 
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Proof. After straightforward calculations, we obtain: 

Bij = — [di{4> o u) ■ dj{ip o n) + dj{(p o u) ■ di^il) o n)) . 

Further calculations shows that in the expansion of 0, as defined in (j3.5p . the third derivatives of 
ip and (j) cancel out: 

= curFcurl [Bij] = dj^B22 + d^^Bn - ^df^Bu 
= d^2{9ii<P°u)-dii'iPou))+df^{d2{^ou)-d2i^ou)) 

— d'l2{di{4> o u) ■ d2{il^ o u) + d2{(t> ou) ■ di{'ip o u)) 
= -{dn{(t>ou) ■ 822^^ o u) + d22{(t> o u) ■ dfiiijou) - 2dl2{4>ou) ■ df2{i^ o u)). 
As a consequence, if S is of class C'^^^'^, the solution of (|3.7p from Theorem 13. II satisfies: 

||i/;n||cfe-i,i < C||0||cfc-i,i < C||(/>||cfc+i,i llV'llcfe+i.i • 
Reasoning as in the proof of Theorem 13. H we obtain: 

||w^tan||cfe.i + ll'"^^llc'=-i>i ^ C'l|-S|lc'=-i>i + C'll</'llc*+i.i IIV'llcfe+i>i ^ C'll'^llcfe+i>i IIV'llc'=+i>i > 

proving the claim. ■ 

4. Spaces of W'^^'^ infinitesimal isometries on developable surfaces 

In this section, we establish some properties of W'^'^ first order infinitesimal isometries on 
developable surfaces S of C^'^ regularity. We give a classification of these displacements and prove 
that they are necessarily C^'^^^ regular. Define: 

V = {V £ W^'^{S,R^); symVV = 0} 

Note that in view of Lemma [2. 5 1 we may freely determine the regularity of any mapping on S, up 
to C^'^ regularity, by considering the regularity of its composition with the chart u. Here and in 
what follows we write / e W'^''^{S) precisely if / o n G Ty^'^(il). 
The following is the main result of this section: 

Theorem 4.1. Let V £V and assume that S is developable of class C'^'^ . Then V G C^'^^'^{S,R^). 
More precisely, writing V = Vtan + {yn)fi we have: 

Vtan G C2'^/2(S,M^) and Vn G C^'^/2(5). 
Moreover, V £V if and only if there exist a,b £ Ty^'^((0, T), M) such that 

(4.1) {Vn){ums,t))) = a{t) + sb{t), 

(4.2) symVVtan{nms,t))) = ^^^l±^Kn{t) {V (t) ^ V (t)) for a.e. {s,t) G M,±, 

1 — SK[t) 

and such that the following integrals are finite: 
(4.3) 

T/ ,N f ( ^ K(a'(t)+sb'(t))\^ dsdt 
Ji{a, b)= b'{t) + -^-^ -f^ — < 00, 

Jm^± V 1 - SK,{t) J I- SK{t) 

(4.4) 

Ma,b)=[ ( a"it) + sb"it) - - smm + '^'^'^[""'^'^ "if^'^^ Xjv^^ < 

Jm ,\ 1 - SK{t) J [I- SK[t)y 
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Proof. Since u G C^'^, from (I2.6p we conclude that k„ is continuous up to the boundary and that 
K is continuous on the open set where Kn differs from zero. 

1. Let V gV. From the observations following Lemma l2. II and from (|2.1ip we deduce that $ 
is bilipschitz on every set of the form: 

M's = eM,±:s€{5- s-{t),s+{t) - 6)} 

with 6 > 0. Indeed, <I>^^ may fail to be globally Lipschitz on <1>(M^±) unless (j2.13p is satisfied. Set 
V3 = (Vn) o u and / = V3 o We have V3 G Ty^'^(O). By a similar reasoning as in Lemma 12.41 
we obtain that <I> is a C^'^ diffeomorphism on M'^. Hence, Lemma 12.51 implies that / G VF^'^(M^) 
with: 

(4.5) dsf{s, t) = VVs{x)N{t), dtfis, t) = {l- sK{t))VV3{x)r'{t) where x = s), 
dlfis,t) = {V^V3ix)Nit))Nit), 

(4.6) dlfis,t) = (1 - SAc(t))(vV3(x)iV(t))r'(t) - ^{t)VV3{x)T'{t), 

dlf{s, t) = (1 - sK{t)f (vV3(x)r'(t))r'(t) + vv^{x) - sK)N{t) - sK'v'it)) . 

Moreover, by ([3.70 : 

(4.7) dlj\s, t) = Vt G [0, T] with K„(t) / 0. 

Indeed, 6 = curl^curl [Bij] = in the present case where symVF = 0. 

Let < ?7 < infigjo,T]{'5^(i), s^(t)} so that {—ri,T]) x (0, T) C Mg±. Denote by /* the precise 
representative of / [6j and define: 

a(t) = r(0,i), b{t) = -{nr],t)-a{t)). 

V 

By (j3.7p the definition of b does not depend on rj, and: 

(4.8) f{s,t) = a{t) + sb{t) V(s,t)GM^±. 

2. Since / G Vl^^'^((0,T) x (—77,17)), foi^ almost every pair si,S2 G {—r],ri) the traces f{-,si) 
and /(•,S2) belong to Vr^'^(0, T), by Pubini's theorem. Hence h = 7j4^(/(-, si) - /(•,S2)) G 
VF2,2((o,r)), and a G VF2,2((o,T)) as well. By Sobolev embedding / G Ci'i/2(M,±). 

Since <I> is a C^'^ diffeomorphism, it also follows that Vn, Vtan S C^'^/^(5') which implies |23j 
that Vtan G C^'^/^. Finally, (j421) fohows from Lemma O 

3. We shall now prove that, given the structure (j4.8p . condition Vn G VF^'^(5) is equivalent to 
a,6 G Ty2,2(-Q^2^) satisfying (gSD, (g^l). This will conclude the proof. 

Inserting into (fOI) we obtain, for ah t G [0, T] \ Jq: 

6'(t) = dlf{s,t) = (1 - ./^)(vV3(x)iV(t))r'(t) - -^(a'(t) + sh'{t)), 
a"{t) + sh"{t) = dlf{s,t) = (1 - sKfN^V^{x)V'{t))V'{t) + k(1 - sK)b{t) - -^^{a'{t) + sb'{t)). 

1 — SK 

Solving for V^V3(x) we get: 

(vV3(x)r'(t))r'(t) = ^ (Q^^(t) + s6"(t) - k(1 - sK)b{t) + -^(a'(t) + 56'(t))), 

(VV3(x)iV'(t))r'(t) = :p^(fo'(i) + T^(«'(0 + 
(vV3(x)7V'(t))7V'(t) = 0. 
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Now a change of variables shows that: 

[ |VV3(x)|2dx= / |(VV3)($(s,t))|2(l-SK)dsdt. 

Jn J M ± 

We see that G W^^'^{n) if and only if ^ and ^ hold. ■ 

We finish this section by pointing out a straightforward corollary of the above calculations: 

Proposition 4.2. Let v £ satisfy 

v{u{^{s,t))) = a{t) + sb{t) for a.e. {s,t)€Ms±. 

Then a,b ^ VF^'^(0,T) and there exists a tangent vector field Vtan G VF^'^(S', M^) to S such that 
Vtan + vneV. 

5. Matching and density of infinitesimal isometries 

Definition 5.1. A one parameter family {us}e>Q C C°'^(5,M^) is said to be a (generalized) Nth 
order infinitesimal isometry if the change of metric induced by is of order e^^^ , that is: 

(5.1) ||(Vn,)^(Vn,) - IdU^is) = 0{e''+') ase^Q. 

Here and in what follows we use the Landau symbols 0{q) and o{q). They denote, respectively, 
any quantity whose quotient with q is uniformly bounded or converges to as g — )• 0. 
Note that if F G V n C°'\ then = id + eV is a (generalized) first order isometry. 

Theorem 5.2. Let S be a developable surface of class C"^^'^, satisfying \2.1^) . Given V G 
V nC'^^^^'^{S), there exists a family {we}e>o C C^'"'^(5, M^), equibounded in C^'^{S), such that for 
all small e > the family: 

Us = id + eV + e^We 
is a (generalized) Nth order isometry of class C^'^. 

Proof. 1. The result is a consequence of the following claim. Let S be of class C^^^'^ and let u^ 
be an (i — l)th order isometry of regularity C^^^'^ of the form: 

j-i 

-u^ = id + ^ e^Wj, Wj G C'^^^'^ 
i=i 

Then there exists Wi G C'^'~^'"'^(S', M'^) so that 0^ = + e^vui is an ith order infinitesimal isometry, 
and: 

i-l 

(5.2) ||wj||cfe-i,i < C*^^ ||t(;j||cfe+i,i 

j=i 

Indeed, setting vui = V £ q2N-i,i g^^^j applying the above result iteratively to find wj G 
Q2N-2j+i,i ^ for j = 2 . . . A^, we obtain the requested = W2 + £w^ + • • • £^~'^wn G C^'^ . 

2. We now prove the claim. Set wq = id. Calculating the change of metric induced by the 
deformation (j)^ we get: 

|(V,/.,)^V</.,-Id 



+ 0{e 



,^+l^ 



MATCHING PROPERTIES OF INFINITESIMAL ISOMETRIES FOR DEVELOPABLE SHELLS 11 

where the expression Aj indicating the change of metric of jth order induced by (pg, is given by: 

j 

Aj = sym((Vwp)^Vwj-p) . 

p=0 

Note that by the assumption Ai = ■ ■ ■ = Ai-i = 0. Consequently, in order for to be an it\i 
order isometry, we must have = or equivalently: 

1 

symVwj = --^ ^ syTaiiVwpfVwi 



2 

p=i 



Applying Theorems 13.11 and 13.21 we obtain that such Wk exists with the estimate: 

i-l 

Wwi^tanWc''^ + ll^j,3|lc'=-i'i — ^ ll^p llc'=+i'i II llc'+i'i ' 

p=l 

provided that all Wp G C^+^^'^ and that S is of class C^"*"^'^. This completes the proof of the claim 
and of the theorem. ■ 

Theorem 5.3. Assume that S is developable, of class C^^^'"*^ up to the boundary, and satisfying 
i2.12\) . Then, for every V €V there exists a sequence Vn G V H C'^'"'^(5, M'^) such that: 

lim \\Vn - V\\w2,2(s) = 0. 

Proof Let a,b G W'^''^{0,T) be as in Proposition Take an,bn G C°°([0,T]) converging in 1^2,2 
to a, b respectively, and define: 

Vn{s,t) = an{t) + Sbn{t). 

By Propositions 14.11 and 14.21 there exist Vn G V such that (V^n) o u o <^ = Vn and \\Vnn — 
Fn||ty2,2(5) — ^ 0. Indeed, the last assertion is equivalent to proving that Jj( 
i = 1,2, which is established immediately after observing that 1 — is bounded away from by 
(I2.13p . Note that (Ki)tan can be chosen suitably such that also || (In)tan — ^ian||vF2.2(5) — )■ 0. In 
view of Lemma 12.41 u G C^'^{Q), that is G C^'^ up to the boundary of S. ■ 

6. The P-limit result 

Consider a family {S^}h>o of thin shells of thickness h around S: 

= {z = p + tn{p); pes, -h/2 < t < h/2}, < h < Hq, 

where Hq is small enough so that the projection map vr : S'^o S,TT{p + tn{p) := p is well-defined. 
For a W'^''^ deformation : ^ M.^ , we assume that its elastic energy (scaled per unit thickness) 
is given by the nonlinear functional: 



h 



5^ 



The stored-energy density function W : W^^^ — > [0, oo] is in an open neighborhood of S0{3), 
and it is assumed to satisfy the conditions of normalization, frame indifference and quadratic 
growth: 

VF G M^^^ Vi? G S0{3) W{R) = 0, W{RF) = W{F), 

W{F) > Cdist\F,S0{3)), 
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with a uniform constant C > 0. The potential W induces the quadratic forms: 

QsiF) = D^W{Id){F,F), Q2{p,Ftan) = uim{Qs{F); {F - F)tan = 0}. 

defined for F G M^^'^, and p G S respectively. Here and in what follows Ftan is the bilinear form 
induced by F on S" through the formula: 

Ftanir, rj) = T ■ F{p)ri £ S, T,ri £ TpS. 

Both forms Q3 and all Q2{p, •) depend only on the symmetric parts of their arguments, with 
respect to which they are positive definite [8j. 

In what follows we shall consider a sequence > such that: 

(6.1) < lim e^//i^ < +00, for some 2 < /3 < 4. 
Also, let: 

/37V = 2 + 2/N. 

We assume that > 1 (the case = 1 is already covered in [23]). Recall the following result: 

Theorem 6.1. [23j. Let S be a surface embedded in R^, which is compact, connected, ori- 
ented, of class C^'^ , and whose boundary dS is the union of finitely many Lipschitz curves. Let 
yh g VF^'^(S''*, M^) be a sequence of deformations whose scaled energies E^{u^)/e^ are uniformly 
bounded. Then there exist a sequence S 50(3) and c^ G such that for the normalized 
rescaled deformations: 

y^{p + tn) = Q^u^{p + h/hotn) - c^ 
defined on the common domain 3^° , the following holds. 

(i) converge in W^''^{S^^'') to tt. 

(ii) The scaled average displacements: 

h /"^o/^ 

(6.2) V^{p) = —= f y^{p + tn) -pdt 

V J -ho/2 

converge (up to a subsequence) in M^^'^(S') to some V £V. 

(iii) liminfft^ol/e''£^^(ti'') >2:(y), where: 

(6-3) T{V) = yJ Q2 {p, (V(^n) - AU)^^^) dp. 

J s 

Here, the matrix field A € W^''^{S,^^'^^) is such that: 

drV{p) = A{p)t and A{p) G so(3) Va.e. peS Vr € TpS. 

In order to prove that the linear bending functional (j6.3p restricted to V is the F-limit of 
the rescaled three dimensional nonlinear elasticity energy {l/e^)E^ we also need to establish the 
limsup counterpart of the F-convergence statement. This is the final contribution of this paper 
which we are formulating in the following theorem. For a full discussion of this topic in this 
context see [9] and [23] . 

Theorem 6.2. Let N > I and assume that S is developable of class C"^^'^ and satisfying \2.12\) . 
Assume that 

(6.4) e^ = o{hl^^). 

Then for every V there exists a sequence E W^''^{S^ ^^M?) such that: 

(i) the rescaled deformations y^{p + tn) = u^{p + th/hon) converge in W^''^{S^'^) to tt. 

(ii) the scaled average displacements given in i6.^) converge in VF^'^(S') to V. 
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(iii) \\mH^ol/e''E^{u'')=X{V). 

Proof. We shall construct a recovery sequence for developable surfaces, based on Theorems 15.31 
and 15.21 Indeed, by the density result and the continuity of the functional X with respect to the 
strong topology of W'^''^{S), we can assume V E VnC^^~^'^(5,M^). In the general case the result 
will then follow through a diagonal argument. 

1. Let e = yf^ jh so e — )• as /i — )• 0, by assumption (j6.ip . Therefore, by Theorem 15.21 there 
exists a sequence w^r : S — > R^, equibounded in C^'^(S'), such that for all small h > 0: 

(6.5) Ue = id + eV + e^w^ 

is a (generalized) A^th order infinitesimal isometry. Note that by (j6.4p we have: 

' - ) - o(/i^+i)//i^+i ^ 0, 



h /iJV+l 



V e 

hence e^"*"-*^ = o{^f?). We may thus replace 0{£^'^'^) with o{'sf^). 

For every p G S, let ^^(p) denote the unit normal vector to Us{S) at the point Usip). Clearly, 
fie e C°'1(5,R3), while by (f53D : 

(6.6) Ue = — 7 = n + eAn + 0{e ). 

Here ri,r2 G TpS* are such that n = ti x T2- Note that since > 1 and is a (generalized) A^th 
order isometry, we have {dnUsl'^ = 1 + 0{e^) and |9ri^£ • dr2Ue\ = C'(^^)) which implies that: 

Using now the Jacobi identity for vector product and the fact that A G so(3), we arrive at (16. 6p . 

Here we introduce the recovery sequence as required by the statement of the theorem. Note 
that the following suggestion for is in accordance with the one used in [7j in the framework of 
the purely nonlinear bending theory for shells, corresponding to the scaling regime (3 = 2. Also, 
a comparison with the similar proof in j 25j for convex shells, with which much of the following 
calculations overlap, is elucidating. Indeed, the main difference here is that instead of an exact 
isometry of the given shell, we make use of an A^th order isometry u^. Consider the sequence of 
deformations £ W^''^{S^,R^) defined by: 

(6.7) u^{p + tn) = Ue{p) + tueip) + -^^d^ip), 

where e depends on h as above. The vector field G Ty^'°°(S', M'^) is taken so that: 

(6.8) lim/i^/^ll 



and: 



(6.9) lim d^" = 2c {p, sym(V(^n) - AU)tan) in L°°(5), 

where c{p, Ftan) denotes the unique vector satisfying Q.2{p, Ftan) = QsiFtan + c®n{p) + n{p) (g) c) 
(see [231 Section 6]). We observe that, as F G C^'^(5, M^) and c depends linearly on its second 
argument, the vector field: 

(6.10) C{p) = c{p, sym(V{An) - AU)tan) 

belongs to L°°(S', M^). Properties (i) and (ii) now easily follow from the uniform bound on Ws and 
the normalization (16.811. 
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2. To prove (hi) it is convenient to perform a change of variables in the energy E^{u^), so to 
express it in terms of the scaled deformation y^. By a straightforward calculation: 

(6.11) -tE\u'') = ^ / f W{VhyHp + tn)) det[Id + th/hoU{p)] dtdp, 

Js J-h„/2 

where Vhy^{p + tn) = Vu^{p + th/hQ-n). We also have: 

Vhy'^ip + tn)n{p) = ^dny'^ip + tn) = n,{p) + th/h^ed^ip), 
h 

(6.12) ^hy^ip + tn)T = Vy^{p + tn) ■ (Id + tn(p))(Id + th / hoU{p))-^ t 

= (Vueip) + th/hoVn,{p) + —h^eVd^p)) (Id + th/hoU{p))-^T, 

for all p G S and r G TpS. 

Prom ([63]), dSS]) and ([ISl) it follows that \\Vhy'^ - Id||ioo(5.ho) as /i ^ 0. By polar 
decomposition theorem, Vhy^ is a product of a proper rotation and the well defined square root 
\/ iy hy^)'^^ hy^ ■ By frame indifference of W we deduce that: 

VF(V^y^) = W (^(WO^VI^^) =w{\d + ^K^ + 0(|i^^p) j , 

where the last equality is obtained by the Taylor expansion, with: 

K"" = {Vhy'')'^Vhy'' - Id. 
As ||-ft^'^||ioo(5'ho) is infinitesimal as /i — )• 0, we can expand W around Id, using the formula: 
1 

W{ld + K) = -D'^W{ld){K,K)+ / {I - s)[D'^W{ld + sK) - D'^W{ld)]{K,K)ds, 
and obtain, in view of using the assumption that W is in a neighborhood of identity: 

(6.13) ^^(V./) = \q, [j^^K^ + -^0{\K^?)^ + ^,o{\K^\\ 

Using (j6.12p we now calculate . We first consider the tangential minor oiK^, as usual conceived 
as a symmetric bilinear form on S: 

KtaniP + tn) = (Id + th/hoU)-^ [id + 0{e^+^) + 2th/ho sym{{Vue)^Vne) 
+ t^h^/hl{Vne)^Vne + o(\/e^)j (Id + th/hoYl)-^ - Id 
= (Id + th/hoU)~'^ 2th/ho sym.{{\/ue)'^\/ne) - 2th/hoU 

+ t^h^/hl{Vne)'^Vne - t^h^/hlU^] (Id + th/hoU)'^ + o{V7'), 

where we used the fact that is a generalized A^th order infinitesimal isometry to see that 
{Vue)'^Vue = Id + 0(e^+i) = Id + o(Ve^), and the identity: 

Ff ^FFf ^ - Id = Ff ^(F - )Ff \ 

By (16. 5p and (16. 6p we also deduce: 

sym((Vue)^Vn£) = H + e sym(V(An) - Ali) + ©(e^), 

(Vn,)^Vn, = n2 + 0(e). 
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Combining these two identities with the expression of K^^^ found above, we conclude that: 
KtaniP + t^) = \/e^(Id + tV^on)"^[2t//iosym(V(^n) -^n)j {U + th / hoU)-^ + o{^/^) . 

Now, as \ne\ = 1, the normal minor of is calculated by means of ()6.12p as: 

n^K^{p + tn)n = |(V/,/)np - 1 = 2th/hoed^ ■ fie + o{Ve^) = 2t/hQ\Q^d^ ■ n + o{\Q''). 
The remaining coefficients of the symmetric matrix K^{p + tn) are, for r G TxS: 

T^K^ip + tn)n = {fie + th/hoed^f(Vue + th/hoVfie + — 2 /i^eVd'') (Id + th/hoUy^r 

\ 



= t/hoVef^idY Vne(Id + th/hoUyW + o(Ve'^), 

where we have used that n^Vn^ = fiJVue = 0. 

3. From the previous computations we finally deduce, with some abuse of notation, that: 

(6.14) iin',^—K'' = ±K{p)tan + ^{C^fi + fi0O inL~(5^°), 

where the vector field C is defined in (j6.10p and the symmetric bilinear form Ktan G L°°(S) is: 

(6.15) K{p)tan = sym(V(An) - AU)tan. 

Using (j6.1ip . (j6.13p . (j6.14p and the dominated convergence theorem, we obtain: 

1 \ r r^o/'^ 

lim -^E^Hv!') = hm^ -h W(Vhy^) det(Id + th/hoIL) dtdp 
h-^o h^o Js 

1 f r^o/s , ^ ^ . 

Q3[-rK{p)tan + T-{C^n + fl(^C)) dtdp 



2. 7 5 7-/40/2 ^^0 
'•ho/2 f2 



1 f r'^°''^ t^ 

- f -^Q2{p,SYm{V{Afi)-AU)tan) dtdp, 

^ Js J-hn/2 "-n 



'S J -ho/2 ri'o 

the last equality following in view of ()6.10p and ()6.15p . Property (iii) now follows, upon integration 
in t in the last integral above. ■ 
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